Quantum Conductance in Semimetallic Bismuth Nanocontacts 
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Electronic transport properties of bismuth nanocontacts are analyzed by means of a low temper- 
ature scanning tunneling microscope. The subquantum steps observed in the conductance versus 
elongation curves give evidence of atomic rearrangements in the contact. The underlying quantum 
nature of the conductance reveals itself through peaks in the conductance histograms. The shape 
of the conductance curves at 77 K is well described by a simple gliding mechanism for the contact 
evolution during elongation. The strikingly different behaviour at 4 K suggests a charge carrier 
transition from light to heavy ones as the contact cross section becomes sufficiently small. 

PACS numbers: 61.16.Ch, 73.50.-h, 73.23.Ad, 73.20.Dx 



In the last few years, several techniques have been ex- 
tensively used to study the conductance of nanometer- 
scale contacts with variable cross sections Most of 
the previous works have been focused on metallic con- 
tacts, for which the Fermi wavelength Xp is very close to 
the dimensions of the atom. The elongation of the con- 
tact produces discrete changes in the cross section (i.e. 
in the number of atoms) resulting in abrupt jumps in the 
experimental G — Z curves (conductance vs. elongation 
of the contact). These jumps are of the order of the con- 
ductance quantum Go = 2e 2 jh since the addition of one 
atom to the contact is accompanied by the opening of 
one or few conduction channels depending on the metal 
p|. This interplay between electronic and mechanical 
properties in metallic contacts makes it difficult to 
observe "pure" quantum effects on the conductance. 

On the other hand, semimetallic nanocontacts could 
provide a unique laboratory to investigate quantum con- 
ductance phenomena in three-dimensional wires. The 
high value of Xf in these systems requires large contact 
areas in order to have total transmissions of the order of a 
conductance quantum. As a consequence, it will be pos- 
sible to distinguish between atomic rearrangements and 
quantum effects. However, earlier experimental results 
did not clarify the situation. The first experiments ||] 
on semimetallic Sb contacts showed several sharp jumps 
in the G — Z curves for conductances much lower than 
the quantum unit, proving that, indeed, many atoms are 
necessary in order to obtain a value of the conductance 
equal to Go- Nevertheless the conductance quantization 
features were absent probably due to a strong backscat- 
tering || or to the large angular opening || of the con- 
tact. More recently, conductance experiments |l0| on 
Bi nanowires at 4K presented a very different behaviour 
showing, in some situations, long conductance plateaus 
close to integer multiples of Go- After statistical aver- 
aging of many consecutive curves, the histogram of con- 
ductance values showed the existence of two peaks close 



to Go and 2 Go- However, the shape of the histogram ex- 
hibited strong deviations from those reported in metallic 
contacts. 
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FIG. 1. Typical G — Z curves corresponding to elonga- 
tion-contraction cycles measured at 77 K: (a) linear scale, 
(b) logarithmic scale. Curves are shifted horizontally for clar- 
ity. 

In order to shine some light on the physics of 
semimetallic contacts, we have performed scanning tun- 
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neling microscopy (STM) experiments on Bi nanocon- 
tacts at different temperatures. First, we analyze in de- 
tail the opening of the first conduction channel, which is a 
particularly difficult task in metallic contacts, where the 
transition from the tunneling regime into contact is gov- 
erned by the well known jump-to-contact phenomenon 
p| . The peculiar electronic structure of bismuth [ p~2|Jl3| 
enables a selection of the carrier type, depending on the 
dimensions of the contact and temperature O] . This ef- 
fect will be studied at 77K and 4.2K for contacts with 
typical conductances in the range < G < 2Gq. The co- 
existence of conductance quantization effects and atomic 
rearrangements will be studied through conductance his- 
tograms of the G — Z curves. Due to the fact that even 
for low conductance values many atoms will be forming 
the contact, the evolution of the G — Z curves will be 
analyzed in terms of a simple gliding mechanism that 
describes the variation of the contact shape during the 
elongation process. 



The experiments have been performed using a high sta- 
bility scanning tunneling microscope (STM) which can 
be operated over the range 2 K - 300 K inside a 4 He 
cryostat. We have used high-purity (99.999%) polycrys- 
tallinc bismuth, cut from a rod, for both STM tip and 
sample. Prior to the measurements the surfaces were 
cut and scratched. At a given location on the sample, 
the tip is crashed forcibly and repeatedly into the sub- 
strate, causing extended plastic deformations that lead 
to a connective neck between tip and sample. Different 
neck geometries can be obtained by varying the indenta- 
tion depth @. In Fig. [l] (Fig. | ) we show typical G - Z 
curves corresponding to a series of consecutive elonga- 
tion cycles measured at 77K ( 4K ). Continuous cycles 
of elongation-contraction of a neck lead to situations in 
which the G — Z curve presents clear steps over the en- 
tire conductance range (typically from to 10 Go). Al- 
though the curves are not reproducible they show a very 
similar shape. Typically, the length of the steps varies 
between 0.2 and 1 nm, depending on the initial condi- 
tions of the elongation-contraction cycle (i.e., the neck) 
|jl5f . These values are close to the lattice parameter of Bi 
(0.475nm). At 77 K the conductance steps have a small 
positive slope without specific structure. However, at 
4K the first plateaus with G > Go usually show smooth 
curved shapes which resemble those observed for metallic 
Al contacts |T^,|l^]. Moreover, the conductance disconti- 
nuities AG are closer to Go, suggesting that at 4K bis- 
muth exhibits a metallic-like behaviour, with a smaller 
Xp. Therefore, only a few atoms are involved in the con- 
tact. Distortions of the lattice parameter at the contact, 
together with the peculiar band structure of bismuth, are 
responsible of the curved plateaus obtained at 4 K. 
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FIG. 2. Typical elongation-contraction cycles measured at 
4 K: (a) and (b) linear scale, (c)logarithmic scale. The cycles 
shown in (a) (the two cycles on the left in (c)) present the 
sub-quantum conductance jumps and steps. The cycle in (b) 
(the rightmost one in (c)) shows a contact breaking process. 
Curves are shifted horizontally for clarity. 



As expected for a semimetallic contact, in the subquan- 
tum conductance regime the G — Z curves show clear 
steps at both temperatures (see Figs. I and 2). Except 
for the small jump discontinuities (caused by rearrange- 
ments of the atoms and thus reflecting mechanical con- 
tact), G decreases exponentially as a function of Z like 
in a vacuum tunneling process, but with an effective tun- 
neling barrier that strongly depends on temperature (see 
Figs, lb and 2b). While at 77K the effective barriers are 
typically of the order of lOmeV, those obtained at 4K 
are close to leV. When the contact finally breaks, the 
conductance is usually so small that, due to limitations 
of our current experimental setup, we are not able to ob- 
serve the transition from "contact" tunneling to actual 
vacuum tunneling. All these facts indicate that, even 
at these very small values of the conductance, tip and 
sample are still in mechanical contact. At 4K there are 
some situations (see Fig. 2(b), and the curve on the right 
hand side in Fig. 2(c)), in which the contact breaks be- 
fore the subquantum conductance regime. In those cases, 
the conductance does not display subquantum steps and 
decreases exponentially as a function elongation, where 
the effective tunneling barrier <j) « 4eV is very close to 
the expected bismuth work function. 
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FIG. 3. Histograms of the conductance values obtained for 
different necks at 77 K. (a) Experimental results from 200 
G — Z curves. The top histogram is an average of histograms 
corresponding to 10 different necks, (b) Histograms obtained 
from the simulation of different initial necks. 



A detailed theoretical description of the electrical prop- 
erties of the contact during elongation is a very difficult 
task. The effective mass of the charge carriers in Bi varies 
over a wide range from m* ~ 0.002too for light electrons 
to to* w too for heavy electrons (mo is the electron mass) 
and depends on the crystal orientation 1 12 |n| . Moreover, 
the peculiar electronic properties of bismuth (scmimetal- 
lic character, band filling, Fermi surface, ...) may also 
vary with lattice deformations due to the extreme uni- 
axial strains occurring in small contacts Jl7t during the 
elongation-compression processes. The apparent metallic 
behaviour at 4K can be qualitatively understood in terms 
of quantum finite size effects |l]|[l4| : as the contact ra- 
dius R decreases the quantized transversal momentum 
k n oc 1/R increases. This leads to a shift in the subband 
energy position E n w h 2 / (2m*)k^ (larger for light 
carriers) and the gap between light and heavy electrons 
increases. At low temperatures, there is a transition from 
light carriers at large cross sections to predominantly 
heavy electron transport for small cross sections where 
the light electron states become empty. This qualitative 
discussion is fully consistent with the analysis of the tun- 
neling regime in mechanical contact (G < Go): The effec- 
tive barriers in mechanical contact are related to the elec- 
tron confinement in the constriction. In a free electron 
picture, the longitudinal momentum of each transversal 
mode or channel n is given by k Zn — \J Kp — (with 
Kp = 2~k/\f) such that to each channel n corresponds 
a ID effective potential barrier in the longitudinal direc- 
tion ||[h|. The modes with k n > Kp — [2n/\p) are 
evanescent or tunneling modes. When R < Xp there is 
no room even for a single channel and the conductance is 
dominated by electron tunneling through the first chan- 
nel with k Zl « ik n w in/R. Assuming that the con- 
striction length increases linearly with the elongation Z 



we will roughly have G oc exp(— j^Z). Therefore, in the 
tunneling regime (before the contact breaks), the slope 
of In G as a function of Z gives the order of magnitude 
of the lateral size of the constriction. Our results then 
suggest, that, at high temperatures, subquantum steps 
are related to structural changes of constrictions involv- 
ing hundreds of atoms (i? is a few nanometers). Since the 
Fermi wavelength is larger or of the order of R, it would 
also be of the order of a few nanometers. In contrast, 
at low temperatures subquantum conductance steps can 
be associated with changes in constrictions consisting of 
a few atoms (R is now a few A). This would imply Xp 
to be one order of magnitude smaller than before. As 
mentioned above, this qualitative discussion is fully con- 
sistent with the behaviour of the G — Z curves in the 
range of a few conductance quanta. 

The observation of subquantum steps in the G — Z 
curves provides evidence of atomic rearrangements in the 
contact. Since these curves are not exactly reproducible 
from cycle to cycle, conductance quantization effects are 
difficult to observe from a single elongation curve. After 
the formation of a given neck, a histogram is obtained 
from 200 consecutive G — Z curves. In fig. ||(a) we 
present a set of histograms obtained from different initial 
necks. As it can be seen, the histograms present peaks 
close to integer multiples of Go which, in contrast with 
metallic contacts, can be unambiguously associated with 
conductance quantization effects. In general, we observe 
that the histogram shape (peak height, number of peaks, 
...) changes from neck to neck p5| . This behaviour is 
clearly different than that observed in some metals (e.g. 
gold nanocontacts [Q). After averaging over 10 different 
necks, the resulting histogram (uppermost in Fig. 3a) 
still shows clear peaks at 1 and 2 conductance quanta 
which now resemble those of metallic contacts. How- 
ever, in semimetallic contacts a variation of a few quanta 
in the conductance requires cross sections much larger 
than those in metallic ones. Therefore, the mechanical 
evolution in these systems is expected to be rather differ- 
ent than in metals. For cross sections larger than a few 
nanometers, yielding occurs as a consequence of collec- 
tive shear events |^,^l| . A qualitative picture of the me- 
chanical evolution can be obtained from a simple model 
system. We shall restrict ourselves to the tensile load- 
ing regime during the elongation of the neck. Suppose 
that the initial constriction is represented by a cylin- 
der of radius Rq and consider the deformation process 
sketched in the inset in Fig. ^(b). The contact is as- 
sumed to glide along a symmetry plane that forms an 
angle 9 with respect to the cylinder axis. The maxi- 
mum elongation Z max before the contact breaks down, 
is given by Z max = 2i?o/tan(#). For a finite elongation 
Z = Z/Z max , the contact area A shrinks according to 
(see inset Fig. £|(b)): 



A(Z) = 2R 2 (acos(i) - Z\/l - Z 2 ^j 



(1) 
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Interestingly, the contact loses its axial symmetry and 
the narrowest section of the contact resembles an ellip- 
tical cross section (Fig. f|(b)) with an axis aspect ratio 
r](Z) w 7T-Rq(1 — Z) 2 /A. An estimation of the trans- 
mission coefficients T nm as a function contact elongation 
can be obtained from a simple model based on the saddle- 
point-contact model [^2|,^3) . The geometry of our contact 
resembles the wide-narrow-wide geometry with 'ellipti- 
cal' cross section discussed in Ref. Q. Therefore, T nm 
can be approximated by: 



the lengths 



1 + exp 

= (n - 
are in 



tana 



(2) 



l/2)/V7 + Vv(m + 1/2) (all 
units of Af)i and tana = 

(A(0) — A(Z)) / ir / Z max defines the effective "opening 
angle" || of the constriction. The conductance curves are 
obtained from the Landauer formula, G = J2nm Tnm- 
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FIG. 4. (a) Experimental G - Z curves at 77 K. (b) Sim- 
ulated G — Z curves according to the model described in the 
text (different initial conditions). Inset: sketch of the defor- 
mation process. 

Typical G — Z curves obtained by this simple model 
for 9 = 50° are shown in Fig. ^ together with some ex- 
perimental results at 77K. In order to mimic the stepped 
structure, the cross section is assumed to be constant 
within a lattice parameter of Bi. Experimentally the con- 
ductance changes from « lOGo to Go after an elongation 
of the contact of the order of 15 nm. Within our model 
this fixes the Fermi wavelength to be in the nanometer 
scale, in agreement with the analysis of the subquantum 



regime. For small Z, i.e. for small opening angles, there 
is almost no scattering (the T„ m 's are or 1). There- 
fore, the conductance is simply given by the number of 
propagating channels in the initial contact No, i.e., those 
nm channels having e nm < \JixA. As the elongation pro- 
ceeds an interesting interplay between the narrowing of 
the constriction and the change of degeneracy of transver- 
sal modes can be observed. For Z = the contact has 
cylindrical symmetry and the transversal levels exhibit 
the well known degeneracy of a parabolic confining po- 
tential. Depending on the initial contact cross section, 
the conductance is then given by 1, 3, 6, 10, . . . Go- As Z 
increases, the cross section decreases and, at the same 
time, the degree of anisotropy r\(Z) also varies, removing 
the mentioned degeneracy. The conductance behaviour 
along the elongation path depends critically on the ini- 
tial conditions of the contact. The theoretical histograms 
(Fig. ||(b)) strongly depend on the initial cross section 
and exhibit an overall shape that is very similar to the 
experimental one. 

In summary, the above experiments provide direct 
experimental evidence for the presence of conductance 
quantization in three-dimensional nanocontacts, as well 
as an experimental exploration of the conductance regime 
prior to the opening of the first quantum channel. Our 
results for moderate-sized contacts are well described by 
a simple model of contact evolution during an elongation 
process that is based on a gliding mechanism. In agree- 
ment with magnetotransport measurements, the differ- 
ences observed between 77K and 4K curves for G < Go 
are consistent with a transition from light to heavy elec- 
trons at the Fermi level when the contact cross section 
becomes sufficiently small. Finally, we consider that our 
work opens a new way to perform experiments on quan- 
tum coherence phenomena in atomic scale contacts that 
up to now were doable with semiconducting heterostruc- 
tures only. 
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